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A black hole window into p -wave dark matter annihilation 
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We present a new method to measure or constrain p-wave-suppressed cross sections for dark 
matter (DM) annihilations inside the steep density spikes induced by supermassive black holes. We 
demonstrate that the high DM densities, together with the increased velocity dispersion, within such 
spikes combine to make thermal p-wave annihilation cross-sections potentially visible in gamma- 
ray observations of the Galactic center (GC). The resulting DM signal is a bright central point 
source with emission originating from DM annihilations in the absence of a detectable spatially- 
extended signal from the halo. We define two simple reference theories of DM with a thermal p- 
wave annihilation cross-section and establish new limits on the combined particle and astrophysical 
parameter space of these models, demonstrating that Fermi is currently sensitive to thermal p-wave 
DM over a wide range of possible scenarios for the DM distribution in the GC. 


Most astrophysical searches for DM annihilation look 
for velocity-independent, or s-wave, cross-sections ( av ). 
Theories with p-wave cross-sections, (av) oc v 2 , have 
largely remained out of reach in standard searches for 
DM. The only previous limits on thermal p-wave DM 
annihilation come from cosmic microwave background 
(CMB) and radio observations, which are sensitive to 
sub-GeV DM annihilating to e + e~ [Tj, while at higher 
masses CMB observations are orders of magnitude away 
from sensitivity to thermal p-wave annihilations [2j. 

Here we show that thermal p-wave DM annihilation 
can be discovered via gamma-ray emission within the 
density spikes that form around the supermassive black 
holes (SMBHs) at the centers of DM halos. The high 
DM densities within such a BH-induced spike, and the in¬ 
creased velocity dispersion required to support such den¬ 
sities, together boost a p-wave-suppressed DM annihila¬ 
tion rate to potentially observable values. Such spikes 
would appear as point sources to gamma-ray telescopes, 
and would contain sharp spectral features strongly indi¬ 
cating a DM origin. As a concrete illustration we consider 
the Milky Way’s SMBH Sgr A*, showing that Fermi’s ob¬ 
servations of the GC already place new constraints on p- 
wave DM annihilation, and opening the door to potential 
discovery of thermal p-wave DM. 

Models of thermal p-wave DM. A thermal p-wave an¬ 
nihilation cross-section is a generic prediction of a broad 
class of well-motivated DM models. If DM is fermionic, 
then in a parity-conserving theory its annihilation to 
spin-0, P-even final states cannot receive any contribu¬ 
tion from the s-wave. Thus, for instance, the well-known 
model of fermionic Higgs portal DM m m exhibits an 
annihilation cross-section that receives its leading contri¬ 
bution in the p-wave. 

We consider here two simple reference models of ther¬ 
mal p-wave dark matter. First is a hidden sector Higgs 
portal (HSHP) model, where a Majorana DM particle 
X annihilates to pairs of dark scalars s. We consider a 
minimal implementation of this model, described by the 
Lagrangian £ = £ kin - \yS (xx + H.c.) + S 2 - S 4 - 


^S' 2 \H\ 2 , where a discrete symmetry S —> —S, x *X 
forbids cubic and linear terms in V ( S ) and a Majorana 
mass term for y. After S obtains a vacuum expectation 
value, it mixes with the Higgs, allowing the mass eigen¬ 
state s to decay to Standard Model (SM) states. The 
four free parameters describing this model can be taken 
to be m x , to *, and the dimensionless Yukawa and portal 
couplings y and e. For m s < 2m^, the branching ratios 
of s are given by the branching ratios of the SM Higgs 
at the same mass, while for m s > 2(2 to*) the decays 
s — i hh ( tt) must also be included. 

Our second example is a hidden sector axion por¬ 
tal (HSAP) model, where Majorana DM annihilates to 
pairs of pseudoscalars a, which can be described by 

the simple Lagrangian £ = £ k i n ~ ^ (XX + H.c.) — 
2 

i\ya(xx~ H.c.) — -vr-a 2 ; if an approximate symmetry 
is responsible for protecting m Q , then higher-order poly¬ 
nomial self-interactions of a are suppressed [45] . The 
pseudoscalar can decay via dimension-five couplings to 
SM gauge bosons. We will consider the case where its 
coupling to gluons is absent at leading order, which has 
the important consequence that Br(a —► 77) is then 0(1). 
For simplicity, we will take Br(a —> 77) = 1, neglecting 
(e.g.) the generic but model-dependent Br(a —» Z 7) that 
opens up for m a > mz■ With this assumption, there are 
again four free parameters describing this model, which 
can be taken to be m x , m a , y , and e = m a / A, where 1/A 
is the dimensionful coupling of the pseudoscalar to pairs 
of photons. 

The range of gamma-ray signatures exhibited by these 
two theories provides a representative guide to the ob¬ 
servability of a general model of thermal p- wave DM 
annihilating ultimately to visible SM particles. The 
thermally-averaged annihilation cross-section in our two 
reference theories may be written 

(av) = Vl ^ zf(z) + 0(V 1 - zv 4 ), (1) 

7Tm x 

where v is the relative velocity, z = m 2 a /m 2 , and \J 1 — z 
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reflects final state phase space. The function f{z) is 


72 - 160z + 165z 2 - 99z 3 + 37 z 4 - f z 5 
~ 32(4 — z) 2 (2 — z) 4 

for the HSHP model, and 


27 r 6 
32 ^ 


( 2 ) 


fa{z) 


(1 ^) 2 

24(2 - z) 4 


(3) 


for the HSAP model. We numerically solve the Boltz¬ 
mann equation using the full thermally-averaged expres¬ 
sion for the cross-section to determine the value of y 
yielding the observed relic abundance, leaving m x and 
z as the free parameters of each model (we assume for 
simplicity that the SM and the dark sector have the same 
temperature at freezeout). This determines the thermal 
annihilation cross-section uthermai as a function of m x 
and z. The coupling e is irrelevant for astrophysical sig¬ 
natures provided that (i) e <C y and (ii) the mediator 
(a, s) decays before BBN. 

These models are both examples of secluded WIMPs 
[5|; their direct detection and collider signals are oc e 2 
and can be parametrically small. The indirect detection 
signal, on the other hand, remains directly tied to the 
relic abundance, and thus for these models gamma rays 
from DM annihilation in a BH spike can easily be the 
most robust route to discovery. 

The Milky Way’s BH spike. A DM density spike 
due to the presence of the central SMBH Sgr A* forms 
inside the radius of gravitational influence of the SMBH, 
r b = M/v q (G = 1). Here M is the mass of the SMBH 
and vq is the velocity dispersion of DM in the halo out¬ 
side the spike. The precise form this spike takes depends 
on both the properties of DM and the formation history 
of the BH, making DM annihilation signals from such 
spikes potentially powerful probes of both the particle 
properties of DM and the evolution of the Milky Way. 

A general DM spike p{r ) may be well-approximated by 
a series of connected power-law profiles. The spike be¬ 
gins its growth from the inner halo, which in a general¬ 
ized NFW halo follows a power law p{r) = p(ro)(r 0 /r) 7c . 
DM-only simulations yield typical values of 0.9 < y c < 
i .2 im ii], while the dissipative collapse of baryons into 
the disk can adiabatically contract the central DM halo 
into a steeper power-law EM!, with values as high as 
7 C ~ 1.6 being possible in the Milky Way El- 

The DM spike grows inside r b « 0 . 277 , (12] , and is 
well-described by p sp (r) = p(r b )(r b /r)' y “ p . The power- 
law index j sp depends on possible formation histories. 
The steepest spikes are formed by the response of a col¬ 
lisionless DM halo to the adiabatic growth of a central 
SMBH, in which case 7 sp = (9 — 27 c )/(4 — 7 C ), which 
for 0 < 7 c < 2 yields 2.25 < < 2.5 |T211Tol . The 

heating of DM from gravitational scattering off of a suf¬ 
ficiently dense and cuspy stellar density within rh could 


substantially soften the DM spike over the lifetime of the 
Milky Way pH] [15] . The final equilibrium spike profile 
attained as a result of this stellar heating has 7 sp = 1.5, 
while a spike that is in the process of being heated would 
have an intermediate value of 7 sp , perhaps 7 sp ~ 1 . 8 ; see 
uni for a recent summary and discussion, including other 
possible spike and halo solutions. 

Once the DM spike attains the “annihilation plateau” 
density p an n = m x /(av)t at r = r.i n , DM annihilations 
become relevant over the lifetime t of the spike (« the 
age of the SMBH). For r < Vi n , annihilations weaken 
the power law growth to p in (r) = p a nn{r in /r ) 1/2 T7] . 
In this innermost region only particles in eccentric orbits 
with apocenters outside contribute significantly to the 
density inside r ln . This weakened profile arises whenever 
the timescale for annihilation in the Galaxy lifetime de¬ 
creases with decreasing r in a canonical spike, which is 
the case for p -wave as well as s-wave annihilations. Fi¬ 
nally, the inner boundary of the spike is at 4 M [18] ,lj9|. 
The resulting density profile thus may be written as [16] 


p(r) =0, r < AM (capture region), 

= 4M < r < n (spike), 

PspH) + Pin{t,r) 

= Pb{r b /r) lc , r b < r < R H (cusp), 

= Ph{Rh/ r ) lH , Rh < r (outer halo), 


(4) 


The dominant contribution to any annihilation signal 
comes from the region where r ~ ?y n . 

Critically, the velocity dispersion increases inside the 
spike to support the power-law increase in density. For 
the (isotropic) velocity dispersion profile, we match an 
approximate, piece-wise continuous solution of the Jeans 
equation in the spike onto a constant dispersion in the 
inner halo, 

\ + J^( 7m ~ 7sp V 

D 11 \ 1 T 7 s P ). 

AM < r < n n (inner spike), (5) 

- — 7 —-, r in <r<— - (outer spike), 

r l + 7s P 1 + 7sp 

= Vq = const, j-p— < r (cusp & outer halo). 

Following m, we adopt the following parameter val¬ 
ues for the Milky Way’s DM halo and SMBH: M = 
4 x 1 O 6 M 0 [201 El], po = 0.3 ± 0.1 GeV cm " 3 [ 22 ], 
v 0 = 105 ± 20 km s " 1 [23], R 0 = 8.46^38 kpc [24], 
and t ann = 10 10 yrs. With these parameters we find a 
spike radius of 77, = 1.7 pc, subtending 0.012°, well be¬ 
low the resolution of current and future gamma-ray tele¬ 
scopes [25ll33j . The remaining parameters of the spike 
solution are the exponents j c , 7 sp ; with no direct mea¬ 
surements of these quantities, we treat them as free pa¬ 
rameters of our model. Typical values for 7 y n fall in the 
range 10 - 3 - 10~ 5 pc. 


J (r) = 


M 1 

r 1 + 7» 
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FIG. 1: Spectra of three possible candidates to contain emis¬ 
sion from a BH-induced DM density spike |34l I35j . In black 
are two example predictions of thermal HSHP DM with dif¬ 
ferent particle masses (determining the shape) and spike and 
halo parameters (controlling the normalization). The solid 
line shows predictions for m x = 45 GeV, m 3 = 10 GeV, 
7 C = 1.3, and 7 sp = 1.8; the dashed line, m x = 110 GeV, 
m s = 50 GeV, and an adiabatic spike in a halo with 7 C = 1.1. 


Observability of p-wave DM: continuum emis¬ 
sion. Our aim here is to demonstrate that there is a 
sizable range of possible spike and halo parameters for 
which emission from thermal HSHP DM annihilation in 
a BH spike at the GC is comparable or greater in bright¬ 
ness to detected gamma-ray point sources in the same re¬ 
gion. Fermi' s Third Point Source Catalogue (3FGL) |34j 
contains several point sources near the GC. The Fermi- 
LATlocalization accuracy for a point source depends on 
its brightness and that of the surrounding diffuse emis¬ 
sion, and is ~ 9 arcmin for the sources in question. The 
Fermi team has associated source 3FGL J1745.6-2859c 
with Sgr A*. This has an integrated flux of 2.18 x 
10~ 8 photons/cm 2 s in the energy range 1-100 GeV. The 
source 3FGL J1745.3-2903c is slightly brighter, with an 
integrated flux $ = 3.87x 10 -8 photons/cm 2 s in the same 
energy range, but is 5.1 arcsec offset from Sgr A*, and so 
is a less likely association. Spectra of these two sources 
appear in Fig. |T] together with spectra from HSHP DM 
annihilating inside possible BH spikes. The remaining 
bright 3FGL source within 30 arcmin of Sgr A* is as¬ 
sociated with a pulsar wind nebula and thus is not a 
BH spike candidate. Fermi's Second Point Source Cat¬ 
alogue [Sr reported a single central point source 2FGL 
J1745.6-2858, previously identified as a BH spike can¬ 
didate m > and shown for comparison. 

In Figs. 2 and 3 we compare the flux from HSHP 
DM to the point source fluxes detected by Fermi. For 
each of the three point sources, we find the minimum 
value of the cross-section such that the primary photon 
flux from DM annihilations exceeds the observed flux in 
any energy bin at more than 95% CL, treating each bin 
as an independent Poisson-distributed variable. We use 



m x (GeV) 



FIG. 2: Sensitivity to HSHP DM from Fermi observations 
of the GC. Top: the minimum value of the DM annihilation 
cross-section for which the primary photon emission exceeds 
observed point source spectra at 95% CL in at least one bin, as 
a function of DM mass, for adiabatic spikes with y c = 1.1 and 
z = 0.2. Bottom: the predicted flux in a box-shaped spectral 
feature from HSHP DM annihilating inside adiabatic spikes 
for 7 C = 1.3 (blue, dashed), 1.2 (cyan, dotted), 1.1 (purple, 
dash-dotted). Results are shown for 2 = 0.99, where the box 
is sufficiently narrow to appear as a line, and compared to 
Fermi's limit [39] on line flux from the GC (black). 


Pythia 8 [35] to generate photon spectra. Scalar branch¬ 
ing ratios in the range m s < 75 GeV are calculated using 
HDECAY for m s > 75 GeV, we use branching ra¬ 
tios from the LHC Higgs Cross-Section Working Group 
[351 . Our results are shown in the top panel of Fig. [2] for 
the representative choice of mass-squared ratio 2 = 0.2. 
The large differences between the maximum cross-section 
allowed by the three different point sources arise be¬ 
cause the flux from a BH spike is only weakly depen¬ 
dent on the annihilation cross-section; for a p- wave spike, 
$ oc (o-/cr t hermai)^ 3_Te ' ! ’^ /, ^ 1+7 ” p ^- This exponent is ~ 1/5 
for adiabatic spikes, and ~ 1/2 for heated spikes. The 
BH spike would outshine detected point sources across 
a wide range of possible spike and halo scenarios, no¬ 
tably including adiabatic spikes for 7 C > 1, as we show 
in Fig. [| 

Greater sensitivity could be obtained in a dedicated 
search using more sophisticated signal and background 
modeling. Signal will also include sizeable secondary 
emission arising from the interaction of DM annihila¬ 
tion products with ambient dust, starlight, and mag¬ 
netic fields. Meanwhile, a given BH spike candidate will 
generically contain astrophysical emission in addition to 
any DM signal. The 3FGL ‘variability index’ for 3FGL 
J1745.6-2859c suggests that it may be time-variable, 
which, if confirmed, would set a floor for the astrophysi- 
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FIG. 3: The region of combined halo+spike parameter space 
where Fermi observations constrain thermal p-wave annihi¬ 
lation cross-sections for m x = 110 GeV. In cyan are con¬ 
tinuum constraints for HSHP DM (z = 0.2), from compari¬ 
son with 2FGL J1745.6-2858, 3FGL J1745.6-2959c, and 3FGL 
J1745.3-2903.c in the dashed, dotted, and dot-dashed lines re¬ 
spectively; the blue solid line shows the line search limit of 
Ref. [39] on HSHP DM (s = 0.99); the purple solid line shows 
the line search limit on HSAP DM (z = 0.99). Dotted blue 
and purple lines show an estimate of the improved sensitivity 
to HSHP and HSAP DM respectively offered by a dedicated 
line search. Regions above the curves exceed observations. 
Adiabatic spikes are indicated by the dotted black line. 


cal contribution to the gamma-ray flux from this source. 

Limits on p-wave DM: box and line searches. 
Given the large systematic uncertainties on DM halo and 
spike distributions, it is difficult to conclusively discover 
or exclude p -wave DM using a continuum signal alone. 
Even for DM with an s-wave annihilation cross-section, 
where support for a DM interpretation of a potential 
gamma-ray signal may be obtained from its extended 
spatial distribution, the subdominant but sharp gamma- 
ray line at E 1 = mr>M remains a smoking gun for a 
DM origin, in contrast to the broad continuum signature 
which may more easily be mimicked by astrophysical pro¬ 
cesses [40] 01]. The analog of a gamma-ray line in our 
reference models is a gamma-ray box , from the decay of a 
(boosted) mediator to a pair of photons 02] . The upper 
and lower endpoints of the box depend on the mass split¬ 
ting between the DM and the mediator, and are given 
by (l ± \/l — z) . In the near-degenerate limit 

z ss 1, the box becomes narrower than Fermi's energy 
resolution, and limits from line searches may be directly 
applied. Sensitivity to the flux in wider boxes in less 
degenerate spectra is ~ 2-5 times weaker than the sensi¬ 
tivity to the flux in line-like features at the same value 
of m x [42103]. 

The line search with the best sensitivity to p-wave 
DM is region R3 from 02], which considers the inner 
3° around the GC and notably does not mask point 


sources. We reinterpret this search as a constraint on 
narrow boxes (z = 0.99) originating from a BH spike 
[46] . We show the resulting exclusions for HSHP DM 
annihilating in adiabatic spikes in the bottom panel of 
Fig. § In this model, Br(s —> 77 ) <10 3 , making the 
current line search less sensitive than the continuum con¬ 
straints; this conclusion would also apply to line searches 
in fermionic Higgs portal models, and to HSAP mod¬ 
els where the pseudoscalar dominantly decays to gluons. 
Thus a DM origin for a potential Milky Way signal in 
Higgs portal models would be established via the dis¬ 
covery of a sharp spectral feature within the emission 
of a previously-discovered point source. For our HSAP 
model, however, the box is the leading signal, resulting 
in much greater sensitivity. Our limits on HSAP DM 
are shown in Fig. [3] In this model the non-optimized 
line search of Ref. 02] is already sufficiently powerful to 
exclude adiabatic spikes given the Galactic parameters 
adopted here. 

As a rough estimate of the potential improvement of¬ 
fered by a dedicated search for lines near Sgr A* we ap¬ 
proximate the gain in significance as \J B3/B0.3, where 
Ho .3 is the background flux in a search region of radius 
of order the angular resolution for high-energy photons, 
0.3°, and B 3 is the background flux in Fermi's search 
region R3. Examining SOURCE class photons with en¬ 
ergies above 10 GeV gives an estimate for this ratio of 
S 3 /H 0.3 « 10. We show the resulting estimate of poten¬ 
tial sensitivity to thermal p -wave DM annihilation in the 
dotted lines in Fig. [3] 

Summary and Conclusions. High densities in DM 
spikes around SMBHs, together with the enhanced veloc¬ 
ity dispersions required to support them, allow p- wave- 
suppressed DM annihilation cross-sections to yield visible 
signals in current gamma-ray telescopes. Using Fermi’s 
observations of the GC, we placed entirely novel con¬ 
straints on thermal p-wave annihilation cross-sections. 
More precisely, we constrained a sizeable range of com¬ 
bined particle and astrophysical models (in much the 
same spirit as Galactic searches for s-wave annihilation 
that depend on a halo model as well as the final state), 
and established a well-motivated range of particle and 
astrophysical parameter space where DM discovery may 
be uniquely possible via the detection of sharp spectral 
features in a central gamma-ray point source. 

Beyond the Milky Way, most bulge galaxies are ex¬ 
pected to host SMBHs, which will in turn create DM 
density spikes of varying steepness. DM annihilation 
within these spikes yields gamma-ray point sources with 
a common spectrum of primary photons, although the 
secondary emission will depend on the local environment 
of each SMBH. This conclusion is true for s-wave as well 
as p -wave DM; the novelty for p -wave DM is that the 
point sources may be observable even in the absence of a 
detectable halo signature. When the point sources are too 
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dim to be resolved, BH spikes provide a novel mechanism 
for p -wave DM to contribute at potentially nonnegligible 
levels to the extragalactic diffuse background. 
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